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The kinetics of collisionless continuous medium is studied in a bounded region on a curved manifold. We have assumed 
that in statistical equilibrium, the probability distribution density depends only on the total energy. It is shown that in 
this case, all the fundamental relations for a multi-dimensional ideal gas in thermal equilibrium hold true. 
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According to Gibbs, the basic object of statistical mechanics is an ensemble of identical Hamilto- 
nian systems. The systems do not interact with each other and the assembly of them makes essentially 
a collisionless continuous medium. From the viewpoint of kinetics, the Hamiltonian systems with elas- 
tic impacts, i. e. billiards, are especially interesting. These are systems where particles move inertially 
inside a bounded region and bounce elastically against the boundaries of the region. As it is shown 
in Refs. [1, 2], in a billiards, the probability distribution density as a function of time t (this function 
satisfies the classic Liouville equation) necessarily has the weak limit as t — > ±oo. This result justifies 
the Zeroth Law of Thermodynamics in the Gibbs theory. The weak limit is a first integral of the 
I/"") ■ Hamilton equations and depends, in the ergodic case, only on the system's energy. It is noted in 
Ref. [3] that this is very often justified even without any ergodic hypothesis : it is important here to 
keep in mind the function class, to which the probability distribution density function belongs. 

In Ref. [3], we developed the thermodynamics of billiards in the Euclidean space. It turns out to 
be possible to extend these observations to the general case of a curved configurational space. 

Let M n be a compact configurational space of a natural mechanical system with n degrees of 
^ ■ freedom, x = (x\, . . . , x n ) be local coordinates on M, and y = (y±, . . . , y n ) be the conjugate canonical 
momenta. The motion to be considered is inertial. So the Hamiltonian is a positively defined quadratic 
form with respect to the momenta: 

If we denote the matrix of coefficients 1 1 (X% j 1 1 ElS A, then H = (Ay,y)/2. 

Let /(•) be a nonnegative summable function of one variable. By analogy with the Gibbs canonical 
distribution, we introduce the density of the stationary probability distribution in the phase space T = 
= T*M : 



P(x,y) 



J J f{(3H)d n yd n x' 

R n M 



Here, the factor f3 is introduced to non-dimensionalize the argument of /. It is customary to take (5 
= l/kr, where k is the Boltzmann constant, and r is the absolute temperature. 
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The denominator in (2) is referred to in Ref. [3] as the generalized statistic integral. It can be 
easily expressed in terms of r, the only external thermodynamical parameter here being M, Riemann's 
volume of the manifold. To this effect, we perform the linear change of variables y i— ► p: 

y = C(x)p, C T AC = E. 

Thus, 



I f[^pf)(detA) 2 d n xd n p _ 



I" M M R" 

where 

v = J (det A~ 1 )^d n x 

M 

is the volume M with respect to Riemannian metric (1), 

n oo 
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r / ( y J dr = const ' 



where T is the Euler gamma function. 

Now we calculate the average kinetic energy: 



where 



a 



' ( ^ ) 



r - , 
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It is interesting to note that the average internal energy of a collisionless medium does not depend on 
the volume, which correlates with Joule's law for ideal gas. 

As it is shown in Ref. [4], if a Hamiltonian is a homogeneous function with respect to momenta, 
then the quantities calculated using the general routines of statistical mechanics and density (2) satisfy 
the First and the Second Laws of Thermodynamics. Let us calculate, for example, the thermodynamic 
entropy. For this, we should first (according to Ref. [3]) write down the following relation: 

which gives the coefficient x. According to the general theory, this coefficient must be a function of 
the statistical integral F. In the case in question, 

_2a J_ 
X briF' 
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Let $(F) be the antiderivative of >c{F). Then, as shown in Ref. [3], the thermodynamic entropy is 
given by the equation 

p d{3 



Hence, 



s = a + 2a lnF = congt + 2a ( lQV + n lnT \ (3) 
b on on V I 



In the case of the Gibbs canonical distribution ( f(z) = e~ z ), one can easily show, upon integration 
by parts, that 2a = nb. 

Usually, the entropy of ideal gas in a three-dimensional vessel II with volume w is 

iVlnu; + ^-hir + const, (4) 

where iV is the number of gas particles (this expression is sometimes multiplied by the Boltzmann 
constant k, but we do without it). To compare (3) and (4), let us consider the Boltzmann-Gibbs gas 
consisting of N identical small balls, moving in a vessel II. The balls collide elastically with each other 
and with the walls of the vessel. Then, obviously, n = 3N, while the volume v is approximately equal 
to w N (for, in the case of non-interacting balls, the configurational space M of the system is the direct 
product of N copies of IT). Having made these remarks, we see that (3) and (4) become identical up 
to the insignificant constant factor 2a/nb that depends on the type of the function /(■) and on the 
number of degrees of freedom in the system. 

On the other hand, the entropy in statistical mechanics is given by the integral 

plnpd n xd n y. (5) 



r 

In the case of the canonical distribution, this integral coincides with the thermodynamical entropy. 
Of course, for more general distributions of the form (2), this remarkable Gibbs' result is not valid. 
However, the Gibbs entropy (5) looks as follows: 

T +lnF, (6) 



b 



where 

n 

2vr2 



7 



oo 

/ rH ~ lf (t) ln/ (t) dr = const - 



We see that (3) and (6) coincide up to an insignificant additive constant and a somewhat less insignif- 
icant constant positive factor. 

The latter remark is very important for the kinetics of a collisionless medium, especially for the 
validation of the Second Law in the case of irreversible processes. The matter is that (as proved in 
Refs. [1] and [5]) if we replace the density p in (5) with its weak limit, then the Gibbs entropy gets 
a nonnegative increment. If the entropies from (3) and (5) were not so closely related, this general 
result would not allow a natural thermodynamical interpretation. 

In the case of ergodic billiards, we can make not only general conclusions on increase in entropy in 
irreversible processes, but we can also calculate these increments. As a simple example, let us consider 
the case where a collisionless medium is initially enclosed in the portion M_ C M (regions M_ 
and M\M_ are separated with a wall), being in statistical equilibrium. After removal of the wall, 
the medium expands irreversibly, tending to fill the whole region M. During this, its internal energy 
(and, consequently, its temperature) does not change. According to (3), the entropy gets a positive 
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increment, proportional to logarithm of the ratio of volumes M_/M + , where M + = M. This result 
is in a good agreement with the predictions of the phenomenological thermodynamics. We should, 
probably, also mention that ergodicity of the Boltzmann-Gibbs gas for a vessel shaped as a rectangular 
parallelepiped was ascertained by Ya. G. Sinai [6]. 

According to Ref. [3], the thermodynamical variable P, conjugate to the volume v, is given as 



Hence, 



P = _I2I. ( 7 ) 

(5 3v K) 



p= 2_akr_ (8) 



This is the equation of state for the considered system in statistical equilibrium. This equation is 
identical in form with the classical Clapeyron equation. If f{z) = e _z , then 2a = nb, and (8) exactly 
fits the Clapeyron equation for a mole of ideal gas. The physical meaning of the variable P is pressure. 

Let us return to the Boltzmann-Gibbs gas of N small balls in a three-dimensional vessel with 
volume w. Then n = 3N, and we can assume that v = w . Substituting this expression into (8), we 
obtain an equation, which is different from the Clapeyron equation. However, there is no contradiction 
here, for P stands for the pressure of 3Y-dimensional gas. The pressure p in ordinary g as , as a 
thermodynamical quantity, conjugate to the volume w, is given by (7), only $ should first be presented 
as a function of r and w : 

_ = _l_cM> dv_ _ 2a kr 

P f3dw pdv dw 3b w ' 

For the Maxwell distribution (where f(z) = e~ z ), 2a/3b = N, and we obtain the classical ideal gas 
equations : 

E = |iVfc-r, pw = Nkr. (9) 

For non-Maxwellian distributions, the value of 2a/36, surely, differs from N. However, within a 
wide range of distributions, for large N, this value is approximately equal to iV: 

2a(N) , N 

lim y —A- = 1. (10) 



N^oc 3Nb(N) 

For example, this range includes distributions with densities 



f(d)= g (r)e- r2 / 2 , (11) 



where g(r) is an arbitrary non-negative polynomial in r. 
Indeed, 



/ 



rn+a - le - r 2 /2dz = _1 / r „ +Q+le _ r 2 /2d ^ 

n + a 



o o 

Since — t — — ► 1 when a is fixed, this results in the limit relation (10). Recall that functions of 

n + a v ' 

the form (11) are referred to as partial sums of Gram-Charlier series, and are commonly used to 
approximate the distribution densities of arbitrary random variables. In the case in question, such 
an approximation is possible due to the well-known observation (traced as far back as to Boltzmann) 
that in the major portion of a high-dimensional space, any distribution is close to normal (strict 
formulations and discussion can be found, for example, in Ref. [7]). Since N, as a rule, is extremely 
large (of the order of 10 23 ) and not precisely known, we can as well use the classical equations (9) 
instead of E = {a/b)kr and pw = (2a/3b)kr. 
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Thus, we have built a complete (nonequilibrium) theory of ideal gas within the framework of 
Gibbs' general approach, using the concept of weak limits of probability distributions and the result 
concerning the ergodic behaviour of the Boltzmann-Gibbs gas. As distinct from Boltzmann's ap- 
proach, we do not use any additional assumptions (like the condition of statistical independence of 
double collisions). The substantial difference from Boltzmann's approach is that in our theory the gas 
reaches statistical (thermal) equilibrium both as t — ► +00 and as t — > —00, these equilibriums being 
identical. This fact is fully consistent with the invertibility property of the equations of motion. 

Note also that by virtue of (10), the entropy equation (3) yields the classical formula (4) for 
monatomic ideal gas. Besides, the statistical entropy (5) coincides (up to an additive constant) with 
the thermodynamical entropy (3) as t — > 00. 

A supplement. Particle distribution functions 



be the distribution density of the Boltzmann-Gibbs gas, which is a system of N small identical balls 
enclosed in a rectangular box; Xj denotes the coordinates and momenta of the j-th ball. The func- 
tion (12) satisfies the Liouville equation and the initial condition p^(x,0) at t = 0. According to 
Bogolyubov (see Ref. [8]), it is useful to introduce s-particle distribution functions p s (x\, . . . ,x s ,t), 
averaging density (12) over x s+ %, . . . , xpf. The particle distribution functions satisfy the infinite chain 
of "hooked" equations, a so-called BBGKY (Bogolyubov, Born, Green, Kirkwood and Yvon) chain. 
Under some additional assumptions (specifically, that of molecular chaos "in the past"), in the case 
of rarefied Boltzmann-Gibbs gas, one derives the kinetic Boltzmann equation for one-particle dis- 
tribution function p\. These assumptions are not self-evident, do not follow from the principles of 
Gibbs' statistical mechanics, and are to certain extent similar to Boltzmann's assumption of statistical 
independence of the balls' velocities before a double collision. Two important facts follow from the 
Boltzmann equation: 

1) Boltzmann's entropy 



monotonously increases with time, and 
2) as t — ► +00, the distribution p\ tends to the Maxwell distribution. 

However, these conclusions (at least, the former) cannot be directly verified by experiment. The 
matter is that the thermodynamical entropy is introduced only for equilibrium states. The ideas of 
determination of entropy for nonequilibrium states (like those given, for example, in Refs. [9,10]) are of 
methodical nature, proposing to introduce an infinite number of additional internal thermodynamical 
parameters. Using these ideas, the reader can find, for example, the entropy of ideal gas as a function 
of time as the gas is adiabatically expanding into vacuum (Joule's classical experiment). 

We develop a different approach in the nonequilibrium statistical mechanics of the Boltzmann- 
Gibbs gas. It has nothing to do with the analysis of additional assumptions that can be used to 
close Bogolyubov's chain of equations. We evolve Gibbs' classical principles and try to avoid entirely 
additional assumptions of conceptual nature. The crucial idea of our approach is: transition to 
thermodynamical (statistical) equilibrium is equal to replacement of the distribution (12) with its 
weak limit. This idea arises very naturally when one proceeds from microscopic to macroscopic 
description of a dynamical system. The weak limit (as t — > +00 and t — > —00) of density (12) (if it 
exists) coincides with Birkhoff's average 



Let 



Pn{xi, ■ ■ ■ ,XN,t) 



(12) 




p(xi, . . .,x N ). 
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Let <p(xi) be a test function. Then 



tppj^d x\ . . . d xn = / tppid xi 

where 

Pi(xi) = I p N d b x 2 ---d b x N . (13) 



If the initial system with 3iV degrees of freedom is ergodic, then p N is a summable function that 
depends only on the total energy: 

PN = f(^(vl + ... + v 2 N ))/v I fd 3 Vl ...d 3 v N . (14) 

R 3JV 

Here m is the mass of the points, vj is the squared velocity of the j-th. ball, the parameter (3 has the 
dimension of the inverse of energy (it is introduced to non-dimensionilize the argument of /) and v is 
the volume of the 3iV-dimensional configurational space of the system of N balls. The denominator 
in (14) makes the integral of pn over the whole phase space equal to unity. 

Formula (14) shows that every possible position of the iV balls is equiprobable. This fact, noted 
earlier in Ref. [2], means that a homogeneous distribution is established in the state of thermal equi- 
librium. In fact, a similar conclusion follows Boltzmann's theory: the density of gas particles gets 
equalized and at the same time the Maxwell velocity distribution is established. 

Hence, the limit one-particle distribution p ± does not depend on the coordinates, and therefore, 
the averaging in (13) can be replaced with merely averaging over the velocities v 2 , ■ ■ ■ , vn- As a result, 
the following simple expression is obtained: 

; f{±(v?+vi + ...+vi))d^...d?v N 

TTD 3iV — 3 

M») = ~ 71 x (15) 

/ f(^-(vl + vl + ... + v%))d^v 1 ...d^v N 

R 3iV ) 

where x = (3/m, u £ R 3 . 

It turns out that when certain additional constraints (of analytical nature, not statistical) are 
imposed on the function /, the limit one-particle distribution function pi tends to the Maxwell dis- 
tribution as N — > oo. That is, for nearly any initial distribution pn(xi, ■ ■ ■ , %N, 0) (even without 
assuming that p^ is symmetrical relative to x±, . . . ,xn), the balls' velocity distribution in the state 
of thermal equilibrium is, to all practical purpose, normal (if, as usual, iV is sufficiently large). 

It should be underlined that, in such an approach, it is meaningless to speak of the rate of 
convergence of p\ (as a function of time) to the limit distribution p 1 , since the p\ itself does not tend 
anywhere at all. One can only speak of the rate of convergence of the average values of the dynamic 
quantities. If, for example, one takes the characteristic function of certain region inside the vessel as 
a test function if, then it will be just reasonable to consider the rate of equalization of the number of 
balls in this region. 

To derive an expression for the limit distribution (as iV — > oo), we put 3iV = m + 2 and trans- 
form (15) : 

rfl + ffi) Jr^f( ul + U ^ + r2 )dr 



= " 7 x • ( 16 ) 
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Here u = (u\, 1x2,113) and T is the gamma function. The density, actually, depends on the velocity \u\ 
Therefore, its k-th. moment 



/ 



/0l(|7i|)|tt| fe dui dU2 dll3 



is equal to 



00 







2r(i + mjr(*±!) jV^V^J* 



(17) 



When deriving this expression, we used (16) and the elementary properties of the gamma function. 

Assume that the variable x takes all real values; then, it is natural to consider the density p\{x) 
an even function. To simplify the notation, we put x = 1 (or replace the function f(z) with f(xz)). 
Our goal is to show that, as m — > 00, a distribution with the density 

2irx 2 p 1 (x), i£l (18) 

tends to the normal one. 

To this end, we accept two assumptions: 

(a) the limit (as m — > 00) density (18) has a finite positive variance (the second moment, k = 2), 
and 

(b) the function / has a summable derivative. 

Condition (a) is the condition of non- degeneracy of the limit distribution, while condition (b) is 
of technical nature and can probably be weakened. Besides, / should decay at infinity faster than any 
power function : otherwise the integrals in (16) and (17) are not defined for every m. 

Putting k = 2 in (17), integrating by parts and using assumption (a), we obtain: 



00 

lim / 2irx 4: pi(x) dx 

m— >oo / 
—00 

00 / C 2 



lim 



J r + 3/ 



f tra+1 t I S 



Jr +1 f (f ) dc 



7e m+3 /' (f ) 

Here 3c is the variance, which exists due to (a). 



7r +3 / (f ) di 

-3 lim = 3c> 0. (19) 
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Let us put / + cf = g. Then, according to (19), 



2 ^ / y v2 

' 

asm-> oo. 

Now, we calculate the limit of the fourth moments (k = 4): 



r +3 s(C) <n / fe'* 3 f(i) <n - o (20) 



lim i ^ . (21) 



o - 2 



According to (20), the integral in the numerator of (21) can be replaced with the integral 



f- 2 \ r f 2 

: m+5 ft I <a_ \ AC — „U _i_ A\ I cm+3 t I S 



r +& /' de = c(m + 4) J r + 7 J df- 

Using (19), it is easy to calculate the limit (21). It is equal to 1 • 3 • 5c 5 . 

In the similar way, we prove that when k = 2n, the limit (17), as m — > oo, is 

(2n + l)!!c n . (22) 

All the odd moments are, obviously, equal to zero. 

Now, let p\ be the density of the normal distribution in the three-dimensional Euclidean space: 

1 



2a 



Then 



1 ~ 
27rx 2 pi(x) = x 2 e 2(7 

V27T(T 3 

Let us calculate the variance of this distribution: 



4 

x -e 2ff dx = 3a 2 . 



/2vra d 

— oo 

The fourth moment is equal to 1 • 3 • 5<r 4 ; more generally, the 2n-th moment is equal to 

(2n + l)!!cr 2n . (23) 

Formulas (22) and (23) coincide if we put c = a 2 . Hence, according to the Chebyshev-Markov 
moment theorem (sec Rcf. [11]), 

lim pi{u) = p(u), «eR 3 ; a = ^fc (24) 
m— »oo 

Up to now, we have been using the assumption that the function / does not depend on the number 
of particles. In general, of course, this is not the case, and instead of a single function /, we have a 
sequence of functions, f m +2- Nevertheless, we can again put g m +2 = fm+2 + c f'm+2 (provided that the 
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limit distribution has a finite positive variance). It is easy to show that (24) remains valid if the limit 
expression (20) is replaced with a more general one: 

OO / oo 

/ r +2fc+ w(f ) di n r +2fc+1 / m+2 (f ) # - o 

' 

as m — > oo for any integer k ^ 1. 

Recall that, for a fixed value of the total energy, the distribution (14) is called microcanonical . 
According to Maxwell and Borel, it can be transformed to the canonical Gibbs distribution if one 
assumes that the total energy is equal to NE, the average energy E of a single particle being indepen- 
dent of N (see Ref. [12]). Of course, this is an important assumption. A more general implementation 
of this idea, applied to an ensemble of weakly interacting identical subsystems, can be found, for 
example, in Ref. [13] . In essence, the author specifies the conditions, under which the microcanonical 
distribution weakly converges to the canonical distribution as the number of subsystems increases in- 
definitely. Besides, as test functions, the author uses so called adders, symmetrical functions of some 
specially chosen canonical variables. Our construction of the normal distribution is based on different 
ideas. 

One should bear in mind that Boltzmann's and Bogolyubov's theories are not free of all these 
problems, either. Suppose that, at the initial time t = 0, the distribution coincides with the 
distribution (14). The distribution (14) is stationary, and corresponds to the state of thermodynamical 
equilibrium (in Gibbs' approach). In Boltzmann's theory, however, density ~p~i(u) (which is given 
by (15)) corresponds, in the general case, to the initial nonstationary distribution and should tend, 
in the course of time, to the Maxwell distribution. In Bogolyubov's theory, we have a similar case: 
not every summable function can be readily used as the density of the initial distribution. It was 
supposed that the velocities in a particle system should in some remote past (when the particles were 
far from each other) be independent (so that every s-particle distribution function was reduced to a 
product of one-particle functions). It should be underlined that this remote past cannot be replaced 
with the remote future (for discussion, see [8]). However, in our approach, the tendency to statistical 
equilibrium is invariant under the time reversal. 

The work was supported by the Russian Foundation for Basic Research (grant 01-01-22004) and 
the Foundation for Leading Scientific Schools (grant 136.2003.1). 
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